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SOME  STUDIES  OF  WHISTLER  INSTABILITY 

by 

F.  W.  Crawford,  J.  C.  Lee  and  J.  A.  Tataronis 

Institute  for  Plasma  Research 
Stanford  University 
Stanford,  California 

ABSTRACT 

This  paper  treats  the  instability  excited  when  a  transverse  electro¬ 
magnetic  plasma  wave,  propagating  in  the  right-hand  polarized  (whistler) 
mode  parallel  to  a  uniform  magnetic  field,  interacts  with  a  stream  of 
monoenergetic  gyrating  electrons.  The  objectives  of  the  work  are  two¬ 
fold,  first  to  obtain  numerical  estimates  of  the  growth  rates  for  various 
parameters,  and  second  to  identify  the  instabilities  as  being  either  ab¬ 
solute  or  convective.  Th«  results  show  that  if  the  background  plasma 
temperature  is  zero,  absolute  instability,  i.e.  temporal  growth,  always 
occurs.  This  instability  may  be  quenched  either  by  collisions  or  nonzero 
electron  temperature  in  the  background  plasma,  leaving  a  convective  insta¬ 
bility,  i.e.  spatial  amplification.  The  absolute  instability  can  also 
be  quenched  by  introducing  a  longitudinal  energy  spread  in  the  beam. 

Even  Wen  whistler  growth  is  predicted  theoretically,  it  may  not  be 
observable  in  practice,  due  to  the  presence  of  competing  instabilities. 

Two  of  these  are  treated  in  the  paper:  perpendicularly  propagating 
cyclotron  harmonic  waves,  and  longitudinal  beam/plasma  interaction. 
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i.  iNTKOttxrrios 


Whistler  Instability  1*  of  interest  in  sress  o  wi*l*l>  *epsr*ted  ** 

magnetosphorlc  physic*,  thermonuclear  fusion,  and  mcroum  beam  plawn* 

amplification.  There  aro  consequently  many  description*  of  It  to  be 

found  in  the  literature,  written  from  differing  point*  of  Mew,  and 

covering  numerous  special  charged  particle  velocity  dial rlbut Ion* . * 

In  this  paper,  wo  prosent  results  of  an  extensive  series  of  computation* 

on  variations  of  the  basic  case  of  a  gyrating  electron  beam  Interacting 

with  a  cold  background  plasma  Only  propagation  parallel  to  the  magnetic 

field  will  be  considered.  A  particular  feature  of  the  work,  lacking  from 

most  previous  studies,  is  a  detailed  stability  analysis  to  determine 

whother  the  instabilities  are  such  as  to  give  wave  growth  in  space  (cun- 

1  o 

vective)  or  time  (absolute) .  In  addition,  potential  competing 

instabilities  have  been  examined  so  that  their  growth  rates  can  be  com¬ 
pared  with  that  in  the  whistler  mode  for  the  same  parameters. 

The  plan  of  the  paper  is  as  follows.  In  Section  2,  we  present  the 
dispersion  relation  for  a  Maxwellian  plasma  penetrated  by  an  electron 
beam  with  both  longitudinal  and  transverse  energy  spread.  Numerical 
solutions  for  several  limiting  cases  of  this  dispersion  relation  are 
examined  in  Section  3  for  real  wave-number  and  complex  frequency,  and 
predict  the  occurrence  of  instability.  The  nature  of  this  instability 
is  discussed  in  Section  4,  where  it  is  shown  that  in  some  parameter 
ranges  temporal  wave  growth  occurs,  and  in  others  spatial  growth.  Com¬ 
peting  instabilities  are  treated  in  Section  5,  and  the  paper  concludes 
with  a  brief  discussion  of  the  results  and  their  implications  in 
Section  6. 
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a  standard  expression  derived  by  combining  ■  *u  col  1  Islonl  ess  tioltrnann 
equation  with  Maxwell's  equations.  For  the  whistler,  the  generalised 
dispersion  relation  is  given  by*J'*6 
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In  this  expression,  the  summation  is  taken  over  all  beam  and  plasma 

electron  components;  ujp  is  the  plasma  frequency;  fQ(v^,v^)  *s  the 

steady  state  velocity  distribution,  and  v  and  v  are  the  components 

II  •*■ 

of  the  velocity  vector  parallel  and  perpendicular  to  the  magnetic  field, 
respectively.  Equation  (1)  is  defined  for  tju ^  <  0  .  For  tp^  >  0  ,  its 
analytic  continuation  must  be  used.  Here  subscript  i  denotes  the 
imaginary  component  of  ip(=  lPr  +  io^)  . 

We  now  specialize  Eq.  (1)  to  the  case  of  stationary  plasma  with  an 
isotropic  Maxwellian  electron  velocity  distribution,  fp(v  and  a 

stream  of  gyrating  electrons  with  velocity  distribution,  fb(v^,v^  -  v  ) 
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3.  NUMERICAL  SOLUTIONS  OF  THE  DISPERSION  RELATION 

Whistler  stability  can  be  determined  by  solving  the  dispersion  rela¬ 
tion  for  ai  complex  and  k  real.  In  particular,  wave  growth  occurs  il 
the  imaginary  part,  iu^(k)  <  0  .  In  this  section,  we  shall  examine  some 
special  cases.  The  type  of  instability  concerned  will  be  characterized 
in  Section  <1. 

3 . 1  Zero  Plasma  and  Beam  Tempt  aturcs. 

As  v  ,  v  -*  0  ,  Eq.  (U)  reduces  to  a  form  first  derived  by  Bell 
*-  clB 

and  Uuncman, 


DCuu.k)  a  cV 
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UU-U)  -kv 
c  0|) 


2  2 
k  <  v  > 

_ i _ 

r 

2(uu-m  -kv  )' 
c  U,| 


=  o  .  (8) 


These  authors  analyzed  Eq.  (8)  in  the  weak  beam  limit  f  luh/uu  )  -•  ol 

1  b'  c  J 

with  uu  -  kv  uu  .  This  implies  that  the  gyrating  beam  electrons 
U||  c 

•see'  the  whistler  frequency,  uu  ,  Doppler-shi fted  up  to  cyclotron 

rcsonanc?  at  uuc  ,  and  may  excite  instability.  As  our  numerical  results 

will  show,  unstable  waves  are  excited  with  essentially  constant  growth 

rate,  over  a  very  broad  spectrum  of  frequency. 

The  way  in  which  the  interaction  occurs  is  illustrated  in  Fig.  1. 

Figure  1(a)  shows  the  uncoupled  whistler,  and  the  line  along  which  the 

Doppler-shifted  frequency  (uu-k  v  )  equals  uu  .  Since  uu  <  uu  ,  it 

On  c  c 

should  be  noted  that  v  <  0,  i.e.,  the  beam  is  directed  opposite  to 

H 

the  direction  of  propagation  of  the  wave.  Figure  1(b)  illustrates  the 

coupled  modes  in  the  absence  of  beam  transverse  energy.  As  k  -  oo  the 

upper  branch  approaches  uu  while  the  lower  approaches  the  beam  line. 

c 

Note  particularly  that  nj  is  real  for  all  k  real.  This  implies  that 

no  instability  is  excited,  and  agrees  with  the  conclusion  of  Neufeld 

and  Wright^"  for  these  conditions.  Figure  1(c)  illustrates  the  effect 

of  beam  transverse  energy.  There  exists  a  value  of  k  below  which  uu 

is  real,  and  above  which  uu  is  complex.  As  k  -  oo  ,  uu(=  uu  +  iu) .  ) 

r  1 


6 


tends  to 


(9) 

If  k  is  real,  there  must  be  a  pair  of  conjugate  complex  roots.  The 

branch  with  uu  >  0  corresponds  to  decay,  and  is  of  no  interest. 

For  the  parameters  of  Fig.  1,  the  maximum  growth  rate  is 

(uu./u,  )  sa  0.019  corresponding  to  a  temporal  growth  rate  of  1.1  dB/ 

cyclotron  period.  The  corresponding  m  and  k  are  in  the  vicinity  of 

the  intersection  of  the  two  modes  shown  in  Fig.  1(a).  Physically,  it 

might  be  expected  that  the  most  unstable  waves  would  occur  near  that 

intersection  since  it  is  there  that  the  gyrating  beam  electrons  see  the 

whistler  propagating  in  the  background  plasma  Doppler-shifted  up  to  u)  . 

On  the  other  hand,  the  gain  is  only  slightly  lower  well  away  from  the 

coupling  region,  where  the  background  plasma  plays  no  role.  This  insta- 
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bility  of  the  beam  alone  has  been  studied  in  detail  by  Bers,  et  al .  , 
using  full  relativistic  theory. 

3 , 2  Nonzero  Plasma  Temperature.  Zero  Beam  Temperature. 

For  the  case  of  v.  =  0  ,  Eq.  (4)  reduces  to 
til 
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Figure  2  shows  the  effect  of  nonzero  electron  temperature  on  the  unstable 
mode  of  Fig.  1(c).  Even  for  (v^/c)  =  0.02,  corresponding  to  a  tempera¬ 
ture  of  102  eV,  there  is  virtually  no  effect  on  the  growth  rate.  Further 
computations  show  that  even  for  (vt/c)  =  0.05,  the  maximum  growth  rate 
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0  .  The 


is  reduced  by  less  than  10  percent  below  its  value  for  v  = 
implication  is  that  in  considering  the  instability,  the  plasma  tempera¬ 
ture  can  generally  be  ignored. 

The  physical  reason  for  this  can  be  made  apparent  by  considering  an 
alternative  electron  velocity  distribution  which  happens  to  give  a  sinpie 
algebraic  dispersion  relation.  Consider  the  resonance  distribution, 
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We  then  obtain  an  expression  similar  to  Eq.  (8)  but  with  (to  —  uu  )  in 

14  C 

the  plasma  term  replaced  by  (to  —  uu  -  ika).  It  is  clear  that  thermal 
effects  can  be  ignored  if  Jto  —  I0c  |  »  |  ka  |  ,  i.e.,  if  the  intersection 
point  between  the  beam  line  and  the  whistler  dispersion  curve  is  below 
the  region  where  the  whistler  alone  would  be  strongly  cyclotron  damped. 

It  is  interesting  to  note  that  by  the  nature  of  the  modification 
made  to  the  dispersion  relation  for  the  case  of  a  resonance  distribution, 
there  are  no  solutions  with  uu  and  k  real  simultaneously.  We  may 
easily  demonstrate  that  this  is  also  true  for  the  Maxwellian.  In  Eq.  (10), 
assume  uu  and  k  real,  and  set  the  real  and  imaginary  parts  equal  zero. 
With  use  of  Eq.  (6),  this  yields 


0 


(12) 


For  uu  ,  uu  ,  v  and  k  finite  and  nonzero,  Eq.  (12)  cannot  be  sat- 
P  t 

isfied.  This  is  an  important  result.  It  implies  that  each  mode  of 
propagation  is  composed  of  waves  that  either  all  grow,  or  all  decay,  in 
time.  Hence,  if  for  some  value  of  v  it  can  be  established  that  a 
mode  has  growing  waves,  then  we  can  conclude  that  every  wave  in  that 
mode  must  grow  with  time,  and  that  this  cannot  be  changed  by  varying 
the  temperature  of  the  background  plasma. 


3 -3  Zero  Plasma  Temperature.  Nonzero  Beam  Temperature. 


We  may  consider  two  cases  here  separately.  First  is  that  of  a  spread 

in  the  beam  transverse  energy  alone.  It  is  clear  from  Eq.  (4)  that  the 

details  of  this  are  unimportant,  since  only  <v^>  appears.  Physically, 

this  is  simply  a  statement  that  all  of  the  beam  electrons  see  the  same 

Doppler-shifted  frequency  (w  -  kv  )  ,  which  depends  only  on  v 

°H  0||  ’ 

and  can  give  up  their  transverse  energy  to  the  cyclotron  resonance.  This 

is  not  so  for  the  second  case:  that  of  a  spread  in  v  Under  these 

°ll 

conditions,  and  with  v^  =  0  ,  we  have  as  the  appropriate  limiting  form 
of  Eq.  (4), 


D(uu,k)  = 


2,  2 

c  k 


2  V 

(JU  +  -  + 

U)  -  w 


uu 


0) 


,V2 


kv 


*11 


<  v2  > 

z(?)  +  — z'cn 

4v2 


=  0 


(13) 


Figure  3  shows  the  dispersion  characteristics  of  the  unstable  mode 
aS  Vt||  varies<  It  wiH  be  seen  that  the  growth  rate  is  very  sensi¬ 
tive  to  this  quantity.  The  instability  is  essentially  quenched  when 
<vt||/c^  =  0.0045,  which  corresponds  to  a  temperature  of  only  5.2  eV 
compared  with  the  beam  energy  of  640  eV. 

It  will  be  observed  from  Fig.  3  that  there  is  a  value  of  k  ,  for 

any  given  (vt„/°)  .  above  which  the  wave  is  damped.  This  value  can  be 
Lll 

estimated  as  follows.  Assume  uu  and  k  real,  substitute  from  Eq.  (6) 
in  Eq.  (13),  and  then  separate  real  and  imaginary  parts.  This  yields, 


U) 


=  0 


,22  2 
k  c  -  u)  + 


2 

V 

U)  -  (U 


+  (1 


>)0!b  =  0 


(14) 

where  y  =  (<v^/2v2  )  ,  and  Eq.  (7)  has  been  used.  Equation  (14)  is 

only  cubic  in  k  ,  and  can  easily  be  solved  to  determine  the  critical 

value  at  which  o)J  =  o  . 

i 
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4. 


INSTABILITY  CLASSIFICATION 


In  Section  2,  we  have  been  considering  waves  propagating  as 

exP  i[jj(k)t  -  kz]  with  k  real  and  ju (k )  satisfying  the  dispersion 

relation.  For  <  0  the  plasma  is  unstable  and,  at  first  sight,  no 

steady-state  conditions  appear  possible.  It  was  pointed  out  several 

18 

years  ago  by  Sturrock  that  there  are  actually  two  distinct  types  of 
instability.  This  distinction  becomes  apparent  only  after  a  spectrum 
of  waves  is  superimposed  by  inversion  of  their  Fourier  and  Laplace 
transforms.  This  inversion  can  be  written  in  the  form 


A(z,t)  =  J  J  §“  A (tu , k )  exp  i(ujt  -  kz)  ,  (15) 

c 

where  A(iu,k)  is  the  spectral  amplitude  of  the  corresponding  wave,  and 
C  is  the  Laplace  contour  in  the  lower  half  complex  w-plane,  parallel 
to  the  real  axis  and  below  all  singularities  of  the  integrand. 

As  t  -•  oo  ,  Eq.  (15)  may  describe  one  of  the  two  forms  shown  in 
Fig.  4.  The  pulse  disturbance  may  grow  in  time  at  every  point  in  space, 
or  it  may  propagate  away  while  growing  in  time,  eventually  leaving  the 
plasma  at  any  given  point  quiescent.  The  former  instability  is  termed 
absolute”,  while  the  latter  is  termed  "convective".  The  importance  of 
making  the  distinction  now  becomes  clear.  If  the  instability  is  convec¬ 
tive,  it  should  be  possible  to  excite  propagating  waves  having  real  i 
and  complex  k  ,  i.e.,  growing  in  space  away  from  the  source.  If  the 
instability  is  absolute,  there  will  be  growth  from  noise  until  a  satur¬ 
ation  condition  is  reached  everywhere, for  which  the  small-signal  theory 
breaks  down. 

This  section  will  be  devoted  to  a  study  of  the  whistler  instabilities 

already  introduced  in  Section  II  to  determine  their  nature.  To  do  so, 

we  will  use  the  now  well-established  method  described  by  Dcrfler, ^  and 
14 

by  Briggs,  which  consists  of  mapping  the  lower  half  of  the  complex 
ju-plano  onto  the  complex  k-planc  by  use  of  the  dispersion  relation.  The 
stability  criterion  then  states  that  an  absolute  instability  is  present 
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if  a  saddle-point  is  formed  in  the  k-plane  by  the  merging  of  two  roots 
of  the  dispersion  relation  which  come  from  opposite  sides  of  the  real 
axis.  If,  in  the  limit,  as  the  imaginary  part  of  uu  tends  to  0  -  ,  a 
root  crosses  the  real  k-axis  without  forming  a  saddle-point,  the  insta¬ 
bility  is  said  to  be  convective. 

4 . 1  Zero  Plasma  and  Beam  Temperatures. 

This  procedure  has  been  carried  out  for  Eq.  (8),  and  typical  results 

are  given  in  Fig.  5(a).  It  will  be  seen  that,  for  the  parameters  chosen, 

a  saddle-point  is  located  at  (kc/u^)  =  7.29-10.31.  The  corresponding 

branch-point  in  the  uj-plane  is  located  at  (il)  /tl)  )  =  0.598-10.062.  The 

s  c 

growth  rate  of  the  absolute  instability  is  given  by  -Oh  ).  which  corres- 

s  i 

ponds  to  a  gain  of  3.4  dB/cyclotron  period. 

The  locus  of  the  branch-point  is  indicated  in  Fig.  6  as  a  function 

of  the  beam  velocity  and  its  transverse  speed,  and  in  Fig.  7(a)  as  a 

2  2 

function  of  beam  density  (uu./uu  )  .  Figure  7(b)  shows  the  locus  of  the 

b  c 

saddle-point  in  the  k-plane  as  a  function  of  beam  density.  These  dia¬ 
grams  were  obtained  by  solving  numerically  the  simultaneous  equations 

D (th ,  k )  =  0  ,  *D^’k)  =  0  ,  (16) 


where  D(it,k)  is  as  given  in  Eq.  (8).  Note  that  the  second  part  of 
Eq.  (16)  is  equivalent  to  (<Vu/dk)  =  0  .  Hence  the  roots  of  Eq.  (16) 
are  the  positions  of  the  saddle-points  in  the  k-plane,  and  the  corres¬ 


ponding  branch-points  in  tie  uu-plane. 

The  curves  in  Fig,  6  show  that,  for  a  given  v  ,  the  growth  rate 

0,1  21/2 

of  the  absolute  instability  can  be  lowered  by  decreasing  <"v^>  .  On 

9  1/2 

the  other  hand,  for  a  given  <vp>  ,  the  growth  rate  tends  to  increase 

as  v  -  0  ,  reaches  a  maximum,  and  then  approaches  a  non-zero  limit 
Oil 

for  v  =  0  .  Figure  7  shows  that  there  is  an  optimum  value  of 
0)| 

(w^/uT)  ,  near  0.2,  which  maximizes  ,  but  that  -k^  increases 

with  increasing  Oju^/hj^)  over  the  range  examined. 

We  now  wish  to  investigate  conditions  under  which  the  absolute 
instability  is  quenched,  to  establish  whether  a  convective  instability 
persists.  Two  possible  mechanisms  for  quenching  the  absolute  Instability 


will  be  considered:  collisions,  and  non-zero  plasma  electron  tempera¬ 
ture  . 

4.2  Effect  of  Collisions. 

Electron-neutral  momentum  transfer  collisions  occurring  with  fre¬ 
quency  v  can  be  accounted  for  by  including  a  collision  term  in  the 
Boltzmann  equation  for  the  plasma  electrons.  The  only  modification 
made  to  Eq,  (8)  is  that  (uu  -  in  the  plasma  term  changes  to 

(<JU  -  iv  -  u)  )  . 

c 

The  locus  of  the  branch-point  as  v  and  <v^>1/2  vary  is  shown 
in  Fig.  8(a).  It  will  be  seen  that,  for  the  parameters  chosen,  for 
(VVP>  /c)  =  0.025,  the  branch-point  crosses  the  real  uu  axis  when 
(v/(JUc )  =  0.033,  For  higher  collision  frequencies,  the  absolute  insta¬ 
bility  is  not  present  since  (uu s ) -^  >  0  •  There  is,  however,  a  convec¬ 
tive  instability,  as  the  k-plane  plot  of  Fig.  5(b)  for  this  case  indi¬ 
cates.  For  the  remaining  two  values  of  (<v2>1/2/c)  in  Fig.  8(a),  very 
high  collision  frequencies  would  be  required  to  stabilize  the  absolute 
instability. 

It  is  interesting  to  examine  how  collisions  modify  the  dispersion 
characteristics.  Figure  8(b)  illustrates  this  for  the  parameters  of 
Fig.  5(b).  There  is  a  weak  maximum  in  ^  ,  near  the  coupling  point 
between  the  beam  line  and  the  plasma  mode,  corresponding  to  a  growth 
rate  of  a  few  dB/wavelength  . 

4 . 3  Effect  of  Plasma  Temperature. 

This  may  be  understood  most  simply  by  choosing  for  the  plasma  elec¬ 
tron  velocity  distribution  a  judicious  approximation  to  the  Maxwelllnn. 
We  shall  again  assume  the  resonance  distribution  of  Eq.  (11),  so  that 
the  dispersion  relation  is  given  by  Eq.  (8)  with  the  denominator  of  the 

plasma  term  changed  from  (ju  -  uu  )  to  (uj  -  ika  -  ij;  )  .  The  treatment 

c  c 

is  then  similar  to  that  for  collisions. 

By  making  (a/c)  sufficiently  large,  the  absolute  instability  can 
be  eliminated,  and  a  convective  instability  can  be  uncovered.  For 
example,  for  the  parameters  of  Fig.  8  and  (Cv^-1  2  c)  =  0.025.  the 
critical  thermal  velocity  to  quench  the  absolute  instability  is  plvcn 
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is  given  by  (a/c)  _  0,0045.  The  spatial  growth  rates  are  then  compar¬ 
able  to  those  indicated  in  Fig.  8(b). 

Computations  for  the  Maxwellian  distribution  have  been  made  using 
Eq.  (10).  The  results  are  shown  in  Fig.  9.  They  illustrate  that  tem¬ 
peratures  above  100  eV  are  necessary  to  stabilize  the  absolute  insta¬ 
bility  for  (<vS>  /2/c)  s:  0.025.  When  the  absolute  instability  is 
stabilized,  a  weak  convective  instability  remains. 

4.4  Effect  of  Beam  Temperature. 

It  was  pointed  out  in  Section  3.3  that  only  a  small  longitudinal 

beam  temperature  is  necessary  to  stabilize  the  absolute  instability. 

Figure  10  has  been  computed  from  Eq.  (13)  and  shows  the  branch-point 

locus  for  three  beam  transverse  energies.  The  ratio  f v2  /cv^l 

tl!/  1  J 

required  for  stabilization  is  approximately  0.04  in  all  three  cases. 

For  thermal  energies  in  excess  of  this,  there  is  weak  convective  insta¬ 
bility,  as  in  the  previous  cases. 
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5.  COMPETING  INSTABILITIES 

In  addition  to  the  whistler  instability  which  may  be  excited  by  a 
gyrating  electron  beam,  at  least  two  longitudinal  wave  instabilities 
may  be  excited.  The  first  is  in  the  cyclotron  harmonic  wave  (CHW)  mode, 
and  needs  only  a  nonzero  transverse  component  of  beam  energy.  The  second 
is  a  beam/plasma  interaction  parallel  to  the  magnetic  field,  ard  needs 
only  a  longitudinal  component  of  beam  energy.  In  this  section  we  shall 
consider  the  growth  characteristics  of  these  two  modes. 


5.1  Cyclotron  Harmonic  Waves  (k  =  0) 

II 

These  waves  are  described  by  the  dispersion  relation^ ’ ^ 


1  + 


2  /. 


dv 


Hn<V 

cn-k  v  -ncn 


=  0  (o>±  <  0)  , 


-L  n=-cn 


z*00  /no)  Sf  „  df  \  „  /k  v  \ 

H  (v  )  =  2 n  f  dv  (—  -  +  k  -  )  J2  ( -i-i)  v 

n  ||  J0  i\vi  ||  bv  j  ny  cx>c  J  1 


(1?) 


where  k^  and  k  (  =  k  in  previous  sections)  are  the  components  of 
the  wave  vector  perpendicular  and  parallel  to  the  magnetic  field,  re¬ 
spectively,  and  the  first  summation  in  Eq.  (17)  extends  over  all  charged 
particle  species  present.  For  k^  =  0,  Eq.  (17)  describes  electrostatic 
wave  propagation  with  no  effects  due  to  the  static  magnetic  field.  For 

k  =0,  we  obtain 

II 


1  ~y  y a  w  — — 

Zw  2  Z-  n'  i7  oo-nco 
(JJ  c 


=  0 


a  (k  ) 
n '  r 
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d  1  9f° 
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~  V  ov 
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(IB) 


for  which  numerous  solutions  have  been  obtained.  Here  we  shall  consider 
two  distributions:  The  Maxwellian,  and  the  ring  distribution  given  by 


fo(V\> 


6(vv>  6(v,rvo,)/ 


2nv 


0A 


(Ring)  .  (19) 
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For  the  Maxwellian,  the  dispersion  relation  is 
,  %  V  exp(-X)  In(\)  r^c 

'  2  Z  l  -  0  >  (20) 

%  n=-®  c 


where  X  h.e  beer  wrlttee  for  (k^/toj2,  end  In  is  .  modified  Bes.el 
function  of  the  first  kind.  Figure  11  shows  the  dispersion  chsrecteris- 
tics  described  by  Eq .  (20)  .  For  a  specified  value  of  (oq2/o)2)  ,  propagat¬ 
ing  branches  lie  between  successive  cyclotron  harmonics,  and^here  is  a 
real  co  solution  for  any  real  k^,  i.e.  the  distribution  is  stable. 

More  generally,  if  J  <  0,  for  all  Vj_  >  0,  this  is  always’ 

true.  The  ring  distribution  is  an  example  for  which  thif  is  not  satis¬ 
fied.  It  has  the  dispersion  relation, 


1 


dp  co-njo 

c 


=  0 


kxv0l 


CD 

C 


(21) 


Figure  12  shows  the  numerical  solution  for  u^kj  for  two  values  of 
(cDp/ug.  If  this  parameter  is  relative  small,  each  propagating  branch 
undulates  about  a  harmonic  of  <d..  As  (o>2/a>2)  increases,  two  successive 
branches  can  couple,  leading  to  regions  in  the  dispersion  characteristics 
where  the  0)  is  complex  for  real  k A,  and  instability  is  present.  The 
imaginary  frequency  components  can  become  very  strong  indeed.  For  example, 
in  Fig.  12  (o^/o^)  is  of  order  unity,  and  growth  rates  of  the  order  of 

50  dB/cyclotron  period  are  implied.  This  figure  is  to  be  compared  with 
the  more  modest  growth  rates  of  about  1  dB/cyclotron  period  found  for 
the  whistler  instability. 


To  compare  electrostatic  instabilities  of  perpendicularly  propa¬ 
gating  ChW,  and  the  electromagnetic  instabilities  of  the  whistler  mode, 
we  shall  consider  a  plasma  containing  a  fraction,  a  ,  ®f  the  electrons  in  a 
Maxwellian  distribution,  and  the  remaining  fraction  (l-tt)  in  a  ring 
electron  velocity  distribution.  This  system  has  the  dispersion  relation, 
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(22) 


where  we  note  that  in  this  subsection  only  uu  refers  to  the  total  elec- 

2  2  P 
tron  density  (=  (w  +  UJ.  )  elsewhere. 

P  “ 

Figure  12  shows  the  threshold  conditions  for  instability  in  the 

first  four  frequency  bands  obtained  numerically  from  Eq.  (22).  When 

a  =  0  ,  all  electrons  are  in  the  ring  group,  instability  first  occurs 

2  2 

in  the  passband  2w  <  uu  <  3m  for  (uj  /(D  )  =  6.62.  In  the  passbands 
c  c  p  c 

0<'U)<iD  uu  <'uu<2uj  and  3m  <"  uj  <  4d  ,  the  threshold  condi- 
c  c  c  c  C  2  2 

tions  for  instability  occur  ,  respectively,  at  (ui  /uu  )  =  17.02,  6.81, 

P  c 

and  6.94  .  As  a  increases,  the  threshold  conditions  become  strongly 

dependent  on  the  velocity  ratio  (\r  /v^)  *  Indeed,  in  some  cases,  the 

threshold  condition  may  be  lowered  if  (v  /v^)  larBe> 

Figures  13(b)  and  13(c)  show  that  instability  is  predicted  for  values 

of  (id2/uu2)  that  correspond  to  stability  when  the  electrons  are  exclu:- 
p'  c 

sively  in  the  ring  group.  This  also  has  the  effect  of  increasing  the 
maximum  attainable  growth  rate  of  the  instability  in  a  given  frequency 
band.  This  is  clearly  illustrated  in  Fig.  14  where  the  maximum  value 
of  ((J^/u)  )  is  plotted  as  a  function  of  a  .  It  will  be  seen  that  the 
growth  rate  of  instabilities  associated  with  the  ring  distribution  alone 
can  be  increased  by  adding  to  that  distribution  a  group  of  Maxwellian 
electrons . 

These  results  suggest  that  CHW  instabilities  might  compete  strongly 

with  the  whistler  instability.  Th  >  electrostatic  instability  will  not 

dominate,  however,  when  a  1  ,  corresponding  to  a  weak  beam  passing 

through  a  plasma.  Figure  13  indicates  that  for  a =»  I  ,  it  is  only 

when  v  »  v  that  instability  can  be  excited,  and  then  only  for  a 
t  2  2 

very  narrow  band  of  (uup/(Dc)  .  For  the  parameter  values  used  in  the 
whistler  computations:  (uu^  +  uu?)/(D2  =  26  ,  a  =  0.96  ,  (v  /c)  =  0.02 

P  DC  “ 

and  (v  /c)  =  0.025  ,  we  have  (v  /v  )  =  1.25,  which  is  well  within 
0i  u  x  x 

the  stable  region  of  the  diagrams  shown  in  Fig.  13.  An  additional 
factor  is  that  the  CHW  growth  depends  strongly  on  f(v  )  ,  whereas  the 
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are  finite  and 


It  is  readily  established  that  Z ,  ^  o  if  uu  and 
nonzero.  Hence,  we  conclude  from  Eq .  (26)  that 


k 

II 


0)  -  k  v 
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i- 1  + 

2  2  ,2 

UU  v.  k 
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Equation  (27)  determines  k  at  which  <ju  is  just  real.  We  have 


(28) 


From  Eq.  (26)  we  see  that  if  Jv^J  <  (U)^ju)|  vQ  |  ,  k*  <  0  ,  imply¬ 
ing  that  the  critical  wave-number  is  purely  imaginary.  This  violates 

our  initial  assumption  that  both  k  and  uu  are  real.  Hence,  when 

II 

this  inequality  is  satisfied,  all  waves  with  k  >  0  must  have  m 

II 

complex.  We  know  that  if  v  =  0  ,  then  there  exists  one  mode  of  pro- 

tll  14 

pagation  which  has  growing  waves.  If  v  >  0  ,  but  within  the  limits 

t|l 

specified  by  the  inequality,  every  wave  belonging  to  that  mode  will  have 
a  complex  frequency  with  <r  0  ,  and  will  grow  with  time.  The  situ- 


k  <  0  ,  imply- 


ation  changes  if  |vtJ  >  vCll  f  >  since  now  k^  >  0  ,  as  defined 

by  Eq.  (27).  Hence,  there  is  one  positive  value  of  k  at  which  <u 

II 

is  real. 

When  the  results  obtained  here  are  compared  with  the  instability 

thre  xold  conditions  of  the  whistler  mode  in  the  same  beam/plasma  model, 

there  appears  to  be  almost  an  inverse  relationship  between  the  two 

cases.  For  example,  when  the  beam  velocity  has  a  thermal  spread  and 

the  background  plasma  is  cold,  waves  with  large  k  are  heavily  damped 

by  a  relatively  small  beam  temperature.  Only  waves  with  small  k  can 

II 

grow  in  time,  but  the  growth  rate  is  almost  zero.  For  the  electrostatic 

waves,  it  appears  that  waves  with  small  k  are  the  first  to  be  damped. 

II 

but  only  if  the  thermal  speed  of  the  beam  electrons  exceeds  a  critical 
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number  defined  by  |  |  =(i^/d)^j  |  .  For  smaller  thermal  speeds, 

all  waves  in  the  mode  grow  with  time.  In  the  whistler  instability 

2  2  2  2 

studies,  we  used  (uu  /u)  )  =  25  ,  (ju./uu  )  =  1  ,  and  (v  c)  =  -  0.05  . 

p'  c  b'  c  0|| 

With  these  values  we  obtain  |v  I  =  0.01c,  which  exceeds  the  thermal 

1  C||  1 

speed  of  0.004c  that  was  found  to  be  needed  to  reduce  the  growth  rate 
of  the  whistler  instability  very  significantly. 

Finally,  we  consider  the  case  where  the  background  plasma  is  hot 
and  the  parallel  beam  temperature  is  zero.  The  dispersion  relation  for 
this  case  is 


(29) 

Assuming  m  and  k  real,  and  setting  the  real  and  imaginary  parts  of 
Eq.  (29)  equal  to  zero, yields 


ID. 


1  + 


2.2 
v  k 
t  ii 


1  + 


ID 

21/2  V  k 
t  I! 


0 


(30) 


Since  Z.  ^  0  if  uu  is  non-zero  and  k  is  finite,  Eq.  (29)  cannot 
1  r  II 

be  satisfied.  Hence,  there  are  no  solutions  of  Eq.  (29)  where  both  <v 

and  k  are  real.  The  implication  of  this  result  is  that  the  insta¬ 
ll 

bility  cannot  be  quenched  simply  by  increasing  the  electron  temperature 
in  the  background  plasma.  It  will  be  recalled  from  Section  3.2  that  the 
whistler  instability  also  has  this  property. 

If  the  plasma  temperature  tends  to  zero,  Eq.  (30)  reduces  to  the 
usual  cold  beam/plasma  interaction  expression  which  predicts  an  infinite 
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Provided  the  beam 


1* 

convective  growth  rate  at  the  plasma  frequency, 
velocity  exceeds  the  plasma  thermal  velocity  bv  an  order  of  magnitude, 

Eq.  (30)  will  predict  growth  rates  typically  of  many  tens  of  dB/wa\e~ 

p  p» 

length  for  the  parameters  used  in  our  whistler  calculations.'  ’ 

To  demonstrate  this,  we  replace  Eq .  (29)  by  t^°  approximate  expression 


(iD 
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solutions  for  which  have  been  given  elsewhere.  Equation  (31)  ma>  be 

solved  for  (k  ./k  )  small  to  give 

V  II1  Hr 


c  (od  /to  )  k  c 
r  v  P  c  I!  1 

t  ,  \  ’  U) 

c  'vA  /c)  c 

On 


(u>  /<jj  )  /  a! 

P  c  (  J 
(vo/c>  \"! 


2\V3 


as  the  conditions  for  maximum  growth.  For  our  usual  parameters,  Eq .  (32) 

gives  (k  c/o)  )  ea  100,  (k  c/ai  )  «  15  (v  /v  )^^  The  growth  rate  is 
Hr  c  ||i  c'  0 1|  t' 

enormously  greater  than  for  the  convective  whistler  instabilities  (see 
for  example  Figs.  7(b),  and  8(b)).  It  may  be  argued  that  since  this 
maximum  growth  occurs  for  o>  the  growth  should  actually  be  computed 

for  (u}/ol>c)  <  1  for  a  strict  comparison.  This  is  not  so,  since  the 
beam/plasma  growth  will  certainly  dominate  the  situation,  and  ultimately 
spread  the  beam  velocity.  This  makes  the  validity  of  the  whistler 
computations  somewhat  hypothetical. 
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6.  DISCUSSION 


The  computations  presented  in  Sections  3  and  U  indicate  that  gyrating 
electron  beams  can  excite  weakly  growing  absolute  instabilities  in  the 
whistlor  mode.  These  can  be  quenched  by  high  collision  frequencies  or 
electron  temperature  in  the  background  plasma,  or  by  a  small  longitudinal 
boam  temperature.  A  weak  convective  instability  remains.  Section  5 
indicates  that,  for  the  saroo  parameters,  strong  longitudinal  beam-'plasma 
interaction  can  be  obtained,  but  growth  in  perpendicularly  propagating 
CHW  is  probably  negligible.  Wo  may  now  consider  these  results  in  rela¬ 
tion  to  laboratory  and  space  plasmas.  Pirst,  as  far  as  fusion  study 
machines  with  beam  Injection  aro  concerned,  it  seems  likely  that  che 
longitudinal  wave  instabilities  will  be  most  Important.  It  should  also 
be  remarked  that  a  more  common  instability  in  such  applications  is 
probably  that  appropriate  to  the  loss-cone  electron  velocity  distribution. 
This  distribution  is  formed  through  loss  of  particles  from  the  ends  of 

magnetic  mirror  machines.  It  is  also  subject  to  longitudinal  wave 

PP 

instabilities.  From  the  point  of  view  of  devices,  it  is  strongly 
implied  that  whistler  amplification  is  too  weak  for  use  in  microwave 
beam/plasma  amplifiers  competitive  vith  currently  available  devices. 

Finally,  we  consider  the  magnetosphere,  for  which  whistler  growth 
mechanisms  have  long  been  postulated  to  explain  phenomena  such  as  VLF 
emissions,  and  very  long  multiple-hop  whistler  trains ,^3*24  In  such 
special  regions  as  the  bow  shock,  where  a  stream  of  charged  particles 
impinges  on  the  earth's  magnetic  field,  there  may  well  be  a  gyrating 
beam  type  distribution.  Elsewhere,  the  distribution  is  more  likely  that 
of  a  high  energy  loss-cone  distribution  interacting  with  a  dense  back¬ 
ground  of  effectively  cold  plasma. 3 "  This  specific  distribution  re¬ 
quires  more  intensive  study  than  it  has  received  so  far  to  cover  the  two 
important  points  emphasized  in  our  treatment  here:  determination  of  the 
nature  of  the  instability,  and  comparison  with  competing  instabilities. 

Several  general  points  may  be  made  in  amplification  of  these  last 
comments.  The  first  covers  the  nature  of  the  instability.  To  maintain 
the  amplitude  of  whistlers  in  very  long  whistler  trains,  some  type  of 
convective  instability  suggests  itself,  the  constant  amplitude  observed 
representing  a  saturation  effect.  It  would  be  very  useful,  then, 


-  21 


to  establish  whether  the  loss-cone  distribution  combined  with  a  cold 
background  plasma  leads  to  convective  growth  in  both  directions  along 
the  magnetic  field  lines,  rather  than  to  absolute  instability.  Absolute 
instability  is  more  likely  to  give  rise  to  continuous  VLF  emissions  of 
the  hiss  type. 

The  existence  of  competing  longitudinal  wave  instabilities  raises 
an  interesting  point  concerning  whistler  amplification  observed  from 
the  ground.  The  growth  may  not  be  in  the  whistler  mode  itself,  but 
in  some  other  mode  to  which  the  whistler  couples  at  entry  and  exit  from 
the  active  region.  This  point  does  not  seem  to  have  received  attention 
in  the  whistler  literature,  A  second  point  concerning  competing  insta¬ 
bilities  is  that  our  comparisons  are  of  growth  rates  based  on  small- 
signal  dispersion  relations.  It  is  also  important  to  establish  the 
saturation  levels,  which  can  be  reached  by  the  various  modes,  and  this 
can  only  be  determined  from  nonlinear  theory.  It  is  possible  that  a 
slower  growing  whistler  mode  might  ultimately  reach  a  higher  amplitude 
than  a  rapidly  growing  longitudinal  instability.  For  the  beam-type 
distribution.it  is  very  likely  that  the  beam/plasma  interaction  along 
the  magnetic  field  lines  will  dominate,  since  many  experiments  have 
shown  that  this  type  of  interaction  can  grow  to  amplitudes  sufficient 
to  disperse  the  beam  completely,  whereas  Section  4.3  indicated  that 
only  a  small  beam  velocity  spread  would  quench  the  whistler  instability. 
It  has  been  suggested  in  the  literature  that  conditions  are  accessible 

in  which  the  small-sjgnal  whistler  growth  would  be  the  stronger.  Whether 

i 

this  is  so,  and  whether  the  large-signal  amplitudes  would  continue  this 
trend,  are  still  open  questions. 
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PLASMA 


FIG  EFFECT  OF  NONZERO  PLASMA  TEMPERATURE  ON  THE  UNSTABLE  MODE 

OF  FIG.  l(p)  [(o£*>p  =  25^  (<D^/lDc)  =  1,  (v0„/c)  =  -  0.05 , 

(<  vV/2/c)  =  0.025]  . 
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FIG.  3.  EFFECT  OF  NONZERO  BEAM  TEMPERATURE  ON  THE  UNSTABLE  MODE  OF 

FIG.  1(c)  [(0>W)  =  25,  (“f/O  =  !,  (vn||/c)  =  -  0.05, 

Q-l/oPC  uu  U 

(<  V2  >  /2/c)  =  0.025]  . 


PULSE  AMPLITUDE 


FIG.  7.  STABILITY  ANALYSIS:  LOCI  OF  BRANCH-POINTS  IN  THE  -PLANE 
AND  SADDLE-POINTS  IN  THE  k-PLANE,  WITH  VARYING 

[(fflV)  =  25,  (v0„/c)  =  -  0.05,  (<  v-  >l77c) 


=  0.025], 
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FIG.  8(a).  STABILITY  ANALYSIS:  EFFECT  OF  COLLISIONS 

[(®W)  =  25,  =  1]  . 


32 


(b)  DISPERSION  CHARACTERISTICS  lw^-0.05) 


(kj.Vf/wc) 


FIG  11  DISPERSION  CHARACTERISTICS  FOR  PERPENDICULARLY  PROPAGATING 
CYCLOTRON  HARMONIC  WAVES:  MAXWELLIAN  ELECTRON  VELOCITY 
DISTRIBUTION. 
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FIG.  12. 


DISPERSION  CHARACTERISTICS  FOR  PERPENDICULARLY  PROPAGATING 
CYCLOTRON  HARMONIC  WAVES:  RING  ELECTRON  VELOCITY 


nxcmTDIITTnW 


FIG.  15  (a)  and  (b) 


PERPENDICULARLY 
A  MIXTURE  OF  (c 
VELOCITY  DISTRI 


FIG.  15  (c)  and  (d).  CRITERIA  FOR  THE  ONSET  OF  INSTABILITY  FOR 
PERPENDICULARLY  PROPAGATING  CYCLOTRON  HARMONIC  WAVES  IN 
A  MIXTURE  OF  (a)  MAXWELL IAN7  AND  (l-a)  RING  ELECTRON 
VELOCITY  DISTRIBUTIONS. 
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This  paper  treats  the  Instability  excited  when  a  transverse  elect romagnetlc 
plasma  wave,  propagating  in  the  right-hand  polarised  (whistler)  mode  parallel  to  a 
uniform  magnetic  field,  Interacts  with  a  stream  of  monoenergetlc  gyrating  electrons. 
The  objectives  of  the  work  are  two-fold,  first  to  obtain  numerical  estimates  of  the 
growth  rates  for  various  parameters,  and  second  to  identify  the  Instabilities  as 
being  either  absolute  or  convective.  The-reeults  show  that  If  the  background  plasms 
temperature  Is  zero,  absolute  instability,  i.e.  temporal  growth,  always  occurs. 

This  Instability  may  be  quenched  either  by  Introducing  s  longitudinal  energy  spread 
In  the  beam.  Even  Then  whistler  growth  Is  predicted  theoretically,  it  may  not  be 
observable  In  practice,  due  to  the  presence  of  competing  instabilities.  Two  of 
these  are  treated  in  the  paper:  perpendicularly  propagsting  cyclotron  harmonic 
waves,  and  longitudinal  besm/plssms  Interaction. 
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